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I. INTRODUCTION 

Advances in computational algorithms as well as hard- 
ware speed have enabled a new era of precision calcula- 
tions in the few-nucleon sector. It is now possible to 
compute binding properties, and in some cases other 
properties, of nuclei up through A=10, starting with the 
Schrodinger equation plus phenomenological two- and 
three-body interactions, to a precision that allows for a 
meaningful comparison between theory and experiment. 

Computations are now precise enough that the quan- 
titative role of special relativity becomes relevant. In 
principle, Poincare invariance is an exact symmetry that 
should be satisfied by all numerical calculations, thereby 
permitting comparisons between theory and experiment 
to rest entirely upon the nuclear dynamics. In practice, 

1. Consistent relativistic computation is much 
more numerically intensive, thus motivating 
non-relativistic calculations in practice. 

2. Most estimates of relativistic effects have been 
quantitatively small enough that the non- 
relativistic calculations have satisfactory precision. 

An additional complication is that a non-relativistic 
model does not imply a unique relativistic extension. 
Relativistic invariance requires the invariance of proba- 
bilities in different inertial coordinate systems; this sym- 
metry can be satisfied in a variety of frameworks and 
models that share the same non-relativistic limit. Since 
information is lost in taking the non-relativistic limit, 
there is no unique way to get it back. 
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At the two-body level, calculations in a single iner- 
tial coordinate system are not constrained by Poincare 
invariance; it only ensures that the results in all other in- 
ertial coordinate systems are equivalent. This affects how 
the two-body dynamics is embedded in the three-body 
problem. In addition, there is a unitary group of three- 
body scattering equivalences that preserves Poincare in- 
variance and leaves the two and three-body scattering 
and bound state observables unchanged, at the expense 
making off-shell modifications and modifications to the 
two and three-body interactions. This freedom is not 
constrained by experiment. This means that the contri- 
butions of relativistic effects, off-shell effects, and three- 
body interactions cannot be uniquely separated in the 
three-body problem. 

Therefore, there is no single, unique "relativistic ef- 
fect" in few-nucleon calculations. However, we find that, 
for a variety of dynamical input and some variation in 
model assumptions, there is a consistent pattern in the 
quantitative aspects of relativity in the calculations. 

There are two primary approaches for modeling rel- 
ativistic few-body problems in quantum mechanics. We 
consider here a class of Poincare invariant quantum mod- 
els 0] of few-particle systems. In these models Poincare 
invariance is an exact symmetry that is realized by a uni- 
tary representation of the Poincare group on the few- 
particle Hilbert space. A number of equivalent represen- 
tations of Poincare invariant few-body models are given 
in 0, Q • Poincare invariant quantum mechanics has as 
its starting point a Hilbert space with a fixed number 
of particles (nucleons) and a set of 2-, 3-,. . . body inter- 
actions. These features it shares with its non-relativistic 
counterpart based upon the Schrodinger equation. It dif- 
fers from the latter in the relationship between interac- 
tions in different inertial reference frames. At the level of 
two-body phenomenology, one can make connections to 
parameters fitted to data on the basis of the Schrodinger 
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equation. If this is done in one inertial coordinate system, 
the Poincare invariance can be used to generate interac- 
tions in any other inertial frame. There are several ways 
to do this, as will be discussed below. 

Glockle, Kamada, and collaborators [USB have stud- 
ied aspects of Poincare invariant quantum mechanics in 
two- and three-body problems. They have introduced 
a mapping between the interactions in the relativistic 
two-body mass operator and the interactions in the non- 
relativistic center-of-momcntum two-body Hamiltonian, 
and have solved the three-body problem for specific in- 
teractions 0- Different methods for utilizing realistic 
nucleon-nucleon interactions in relativistic calculations 
were given by Coester, Pieper and Serduke || and by 
Glockle, Lee, and Coester 9]. Our goal here is to identify 
the dominant sources of relativistic effects for a variety 
of two-body mappings and interactions. 

The second approach to the relativistic few-body prob- 
lem uses quasipotential equations. These are relations 
between covariant amplitudes in local field theory. When 
some of the amplitudes are treated as input, these re- 
lations become equations for the remaining amplitudes. 
Matrix elements of few-body observables in eigenstates 
of the four momentum can be calculated from the solu- 
tion of the quasipotential equations using Mandelstam's 
method [jjj, . For systems of strongly interacting par- 
ticles, the input to these equations is not known and must 
be modeled, as in Poincare invariant quantum mechanics. 
While calculations of the same observables can be com- 
puted in both approaches, there exist no unique relation 
between quasipotential equations and Poincare invariant 
quantum mechanics. Nevertheless, quasipotential mod- 
els have played a historically important and valuable role 
in motivating the structure of model nucleon-nucleon in- 
teractions. The nature of the relativistic effects depends 
upon the specific assumptions used to extract quantum 
mechanical interactions from the quasipotential equa- 
tions. Stadler and Gross have studied solutions to the 
three-nucleon bound state problem in the context of the 
spectator approximation to a meson-nucleon field the- 
ory. Sammarruca and Machleidt |13j have studied spe- 
cific kinematic factors that arise in extracting quantum 
mechanical interactions from quasipotential equations in 
an effort to identify the dominant sources of relativistic 
effects. 

In this paper, we attempt to identify the important 
quantitative relativistic effects within the framework of 
Poincare invariant quantum mechanics. The goal in this 
paper is not to provide complete solutions to a set of 
three-body problems, but rather to understand where 
the major relativistic effects occur. In fact, the most 
important effects are embedded in multiplicative kine- 
matic factors, which makes it relatively easy to incor- 
porate them into a Schrodinger-based momentum-space 
three-body code. 



II. THE TWO-BODY PROBLEM 

One-nucleon momentum/spin eigenstates satisfy the 
normalization condition: 



(p',m'IP)M) = S(p' - p)<5 p 



(1) 



These states transform as mass m spin 1/2 irreducible 
representations of the Poincare group 3]. For the two- 
body problem it is useful to use a basis that transforms 
irreducibly with respect to the tensor product of two one- 
body representations. These Poincare irreducible eigen- 
states have the structure 



\P,j,/j,,k;l,s) 



(2) 



where P is the total linear momentum, j is the canonical 
spin, (i is the z-component of the canonical spin, and 
k is related to the invariant invariant mass of the non- 
interacting two-body system by 



M {0) = 2y/m 2 + k 2 . 



(3) 



The quantum numbers I and s are degeneracy quantum 
numbers that determine the multiplicity of each represen- 
tation with the same and j; I and s have the same 
spectrum as the non-relativistic total spin and orbital 
angular momentum quantum numbers. The overlap co- 
efficients with the tensor product of single particle states 
are Poincare Clebsch-Gordan coefficients which can be 
found in ( [13 Q 0| ) - 

The two-body mass operator is defined by adding in- 
teractions to the non-interacting mass operator: 



M = M (0) + v. 



(4) 



In Poincare invariant quantum mechanics the delta func- 
tions that multiply interactions are important. The two- 
body interaction acting on the two-body Hilbert space is 
denoted by v; v denotes the internal two-body interaction 
related to v by 

(P',f, /A k 1 , 1', s'\v\P,j, k, I, s) = 



S(P' -P)S j/j 6 l± ^(k',l\s'\v j \k,l,s). 



(5) 



With this choice of v it is possible to find simultaneous 
eigenstates of M, P, and • z. These states trans- 
form like mass M, spin irreducible representations of 
the Poincare group. Since these eigenstates are complete 
they define the interacting two-body dynamics. 

Two-particle scattering is described by the transition 
operator 



t{e) = v + v 



1 



e -MO — M 



(G) 



which satisfies the relativistic Lippmann-Schwinger equa- 
tion: 



t(e) = v + v 



1 



s + iO - Af(°) 



t(e). 



(7) 
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This transition operator has the same relation to the 
differential cross section as the corresponding non- 
relativistic expression, except that the single particle 
masses, m, appearing in the incident current and phase 
space factors are replaced by the corresponding single 
particle energies, oj — \/k 2 + m 2 . A proof based on time- 
dependent scattering theory is given in Q . 

The reduced two-body transition operator t(e) is re- 
lated to t(e) by 

p',j',li!,k',l'8'\t(e)\P,j,ii,k,l,8) = 



5(P' - PfcV/.^ I', s'\t j (e)\k, I, s). (8) 

III. CONNECTIONS TO TWO-BODY 
PHENOMENOLOGY 

All realistic descriptions of the quantum mechanical 
two-nucleon system utilize adjustable parameters in the 
interactions to fit published nuclcon-nucleon phase-shift 
data and the deuteron binding energy. For Poincare in- 
variant quantum mechanical models, these parameters 
involve interaction strengths and ranges; for quasipoten- 
tial models, the parameters may include meson-nucleon 
coupling strengths and form factors. Since there is an 
extensive investment of effort in fitting nucleon-nucleon 
data to solutions to the Schrodinger equation, it is advan- 
tageous to directly use existing high-quality interactions 
to construct equivalent relativistic interactions. 

The relationship between relativistic and non- 
relativistic models fit to the same data is more com- 
plicated than the relation obtained by taking the non- 
relativistic limit of the relativistic model. Understand- 
ing the nature of the fitting process is the first element 
needed to understand the nature of relativistic correc- 
tions. The scattering cross section is a relativistic in- 
variant [T?) ; however, the angular distributions and spins 
have a frame dependence with known kinematic transfor- 
mation properties. Kinematic Lorentz transformations 
are used to correctly transform measured laboratory dif- 
ferential scattering cross sections to the center of momen- 
tum frame. The properly transformed data are used to 
fit interactions that reproduce this data by solving the 
non-relativistic Schrodinger equation. 

The scattering operator and consequently the phase 
shifts for a given angular momentum are Poincare invari- 
ant. The phase shifts can be tabulated as functions of 
the invariant center of momentum momentum or energy, 
5j(k) or Sj(e). These functions are equal when the center 
of momentum energy e and invariant momentum k are 
related by the correct relativistic relation. For equal mass 
particles this relation is e = 2\Jk 2 + m 2 which implies 

Si (A;) = % (2Vm 2 + k 2 ) 5' 3 (e) = Sj - m^j . 

(9) 



In the fitting procedure these phase shifts are identi- 
fied with phase shifts obtained by solving the Schrodinger 
equation. This can be done in one of two inequivalent 
ways [18J. The first is by identifying Sj(k) with the cor- 
responding fc-dependent Schrodinger phase shift and the 
second is by identifying Sj(e) with the corresponding e- 
dependent Schrodinger phase shift. These are inequiva- 
lent because non-relativistically e and k are related by 
e = k 2 jm + 2m. 

If these procedures are used there are no relativistic 
corrections in the center of momentum frame when com- 
pared to a relativistic model fit to the same data. The 
Schrodinger equation produces the exact Poincare invari- 
ant phase shifts either as functions of the center of mo- 
mentum energy or momentum (but not both). 

The problem that can arise is if the phase shifts are fit 
to transformed cross section data as functions of k and 
the potential is constructed by matching the phase shifts 
obtained from the Schrodinger equation as functions of 
e, then the phase shift or cross section predicted at en- 
ergy e = k 2 /m + 2m will be the same as the "measured" 
phase shift a shifted energy e' = 2\/m 2 + k 2 . Similarly, 
if the phase shifts are fit to transformed cross section 
data as functions of e and the potential is constructed by 
matching the phase shifts obtained from the Schrodinger 
equation as functions of k then the phase shift or cross 
section predicted at momentum k = \J em — 2m 2 will be 
the same as the measu red phase shift a shifted momen- 
tum k' = \J e 2 /4 — m 2 . 

This inconsistency increases with momentum. At rela- 
tivistic energies the inconsistencies are of the same scale 
as typical relativistic corrections. This is illustrated in 
Figure 1, which compares the non-relativistic and rela- 
tivistic energies e nr (k) to e r (k) for a nucleon with mo- 
menta k up to 1 GeV. These curves show the scale of the 
energy mismatch as a function of k ; it is a 7% effect at 
1 GeV. Figure 2 illustrates the same effect with the total 
cross section. The curves in Figure 2 compare the total 
cross sections in the Born approximation for the Malfliet- 
Tjon V potential [lj|- The curve in Figure 2 is the ratio 
a{k)/a{k^l + k 2 /Am 2 ) for k up to 1 GeV/c. In this ex- 
ample figure shows that an inconsistent treatment of the 
phase equivalence, while small at low momenta, leads to 
a 15% percent error in the total cross section at 1 GeV/c. 

Thus, the first question that needs to be considered is 
how the non-relativistic interaction was constructed. For 
Argonne U18, which is a typical example of a realistic 
nucleon-nucleon interaction, the phase shifts are deter- 
mined as functions of laboratory beam energy |2f)j | which 
is converted to center-of-momentum momentum using 

~Pb ■ Pt = e-Labm = (k 2 + m 2 ) + k 2 . (10) 

Thus, for the AU18 interaction, the solution of the 
Schrodinger equation reproduces the invariant phase 
shifts Sj(k) as functions of the invariant momentum. 

The methods discussed in this paper make use of the 
phase equivalence between the interactions in the rela- 
tivistic and non-relativistic Schrodinger equations. We 
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A. Coester-Pieper-Serduke (CPS) 




k(GeV/c) 



This method was formulated to study relativistic ef- 
fects in nuclear matter. It uses a nucleon-nucleon interac- 
tion, such as AV18, fit to Sj(k) to construct a scattering 
equivalent relativistic interaction. 

CPS add an interaction to the square of the non- 
interacting mass operator to get a Schrodinger-like equa- 
tion: 



M 2 = M (0)2 + u = 4mh 



k 2 u 
h:= — + — 

m 4m 



11) 



The operator u/4m — v nr is identified with the inter- 
action that appears in the non-relativistic Schrodinger 
Hamiltonian, h. 

The relationship between u in Eq. and the rela- 
tivistic interaction v in Eq. J^J is 



v 2 + {M ( -°\v}. 



(12) 



FIG. 1: Relativistic vs Non-relativistic CM energy as a func- 
tion of k. 




It would be tedious to evaluate v directly; however, it is 
never needed because M, M 2 and h = k 2 /m + v nr , all 
have the same eigenstates. 

If t(e) is generated by the "Lippmann-Schwinger" 
equation 



t(e) = u + u 



1 



e 2 _ M (0)2 + iQ + 



t(e), 



(13) 



and IV'fe") is a scattering eigenstate with momentum k 
and \ipk) is the corresponding plane wave state, then the 
following half-shell relations follow: 



t(e(k))\^ k )=v\^) = (M-M )\^) 



{M 2 -M 2 )\^ k 



+\ - 



e{k) + M 



k (GeV) 



s(fc) + Mo 



e(k) + M 



t{e(k))\t/> k ) 



e(fe) + Mo e(fc) + Af 

FIG. 2: Differences in invariant cross section identified as 

functions of relative momentum to relative energy. Taking matrix elements in the irreducible plane wave 

states \k,j,l,s) gives the relation 

(k',l', S '\P(e(k))\k,l, S ) = 



;(fc) + e^'))]" 1 ^', I', s'\P(e(k))\k, I, , 



consider three ways of utilizing two-body phenomenology 
based upon the fitted interactions in Schrodinger equa- 
tion in Poincare invariant quantum mechanics. These 
methods relate the half-shell relativistic two-body tran- 
sition matrix elements to the half-shell non-relativistic 

transition matrix elements under different sets of assump- 4m , , , _■ 2 . . . 

tions about the phase shift fitting. e (k) + e(fc') ' ' S 1 nr( /m " ' ' S ' ^ ' 
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where e(fc) = 2ui(k) — 2\Jk 2 + m 2 , e nr {k) = k 2 /m, and 
i? is the reduced non-relativistic two-body transition op- 
erator. 

This gives the desired relation between the relativistic 
and non-relativistic half off-shell transition matrix ele- 
ments. This formula will be used to construct the kernel 
of the three-body equations. 

It has been shown that for the CPS method the rel- 
ativistic and non-relativistic cross sections are identical 
functions of the invariant momentum k and the rel- 
ativistic and non-relativistic bound states have the same 
bound-state wave numbers. This phase equivalence is 
exact. 



B. Glockle-Lee-Coester (GLC) 

Glockle, Lee and Coester introduced an approximate 
variation of the CPS method that did not require intro- 
ducing the interaction u. This approximation was used 
in the first calculation of a three-body binding energy 
based on Poincare invariant quantum mechanics |9J. 

Like the CPS method, the GLC method uses a nucleon- 
nucleon interaction, such as ^4V18 [2(j, fit to 8j(k) to con- 
struct a scattering equivalent relativistic interaction. In 
this case the scattering equivalence is only approximate. 

GLC 9] define a new interaction: 



of (jTSj) by the non- 



(k',l',s'\vi r \k,l,s) 



^ (k',l',s'\v j \k,l,s)*l—, (16) 
m V m 



with the corresponding t matrix elements: 
(k',l',s'\ii r {k 2 /m)\k,l, S ) := 



^ (k',l',s'\P(e(k))\k,l,s)J-. 
m V tti 



(17) 



Equation lfT7|l replaces (fTK|) in the CPS method. 
With the above definitions if i satisfies the relativistic 
Lippmann-Schwinger equation, it does not follow that 
t nr will satisfy the non-relativistic Lippmann Schwinger 
equation, however is straightforward to derive 

(k',l\s'\ti r (k 2 /m)\k,l, S ) = (k / ,l / ,s , \vUk,l,s) + 



E 

l"s" 



k" 2 dk" ' ' s 'l u " r ^"' ' ' 



k 2 /m- k" 2 /m + i0 



(k",l",s"\ti r (k 2 /m)\k,l,s) 



1 



2a/ 



(18) 



The CPS approximation is equivalent to neglecting the 
term, " 2 ~y, , which is zero on shell. Neglecting this term 
gives the non-relativistic Lippmann-Schwinger equation. 



CPS approximate the solution, 
relativistic transition operator. 

It can be shown that this method does not lead to 
an exact phase equivalence, it is possible to construct 
systematic corrections that converge to the CPS result. 
In the original CPS calculation this approximation was 
found to be accurate. It was improved using a slight 
adjustment of the interaction parameters. 



C. Glockle-Kamada (GK) 

The GK approach is designed to produce a relativis- 
tic dynamical model that is phase equivalent to a corre- 
sponding non-relativistic model where the phases shifts 
are fit as functions of center of momentum energy, 8'{E). 

The GK approach uses a unitary rescaling of the 
momentum variables to change the non-relativistic ki- 
netic energy into the relativistic kinetic energy. 

For each value of the momentum k, GK define the 
momentum q by identifying the relativistic and non- 
relativistic energies 



2m + — = 2\/ m 2 + k 2 . 



rn 



(19) 



This can solved for be cither for k(q) or q(k). Defining 
h(q) as 



h 2 (q) 



q 2 dq 
k 2 dk 



(20) 



Glockle and Kamada identify a "relativistic" interaction: 
(k',l',s'\v j \k ) l ) s) := 

h{k'){q{k'),l',s'\vi r \q{k),l^)h{k) (21) 
and t matrix: 

(k',l',s'\P(e(k))\k,l,s) := 

h(k')(q(k')J', s'\ii r (e nr (q(k))\q(k), I, s)h(k). (22) 

Then t 3 and i J nr satisfy the relativistic and non- 
relativistic Lippmann-Schwinger equations, as functions 
of k or q, respectively. 

With this method J52J) replaces JT7|) or (|T5jl . 

In Kamada and Glockle show that with the KG 
method the relativistic and non-relativistic phase shifts 
are identified as functions of the invariant energy E, and 
the relativistic and non-relativistic bound states have 
identical binding energy. 
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IV. THE THREE-BODY PROBLEM 

To construct kinematic variables for the three-body 
system, let denote the invariant mass of the non- 

interacting three-body system and be the invariant 
mass for the non-interacting two-body sub-system con- 
sisting of particles i and j. 

Plane-wave basis states are tensor products of one- 
body Poincare irreducible representation space basis vec- 
tors 



|Pl,2 / l,P2,^2,P3,f3) 



(23) 



The two-body interaction V k i and reduced two-body 
interaction Vki are defined in this basis by 

(P,i,/x, (kl)(m)\V kl \P',j\fj,', (kl)(m)') = 



S(P-P')S jr 5^,((kl)(m)\V kl \(kl)(my) (30) 

where 

((kl)(m)\V kl \{kl)(m)') = S{K "~? kl) x 

K k iK kl 



As in the non-relativistic case it is useful to define the 
total momentum 



i=i 



(24) 



and relativistic Jacobi momenta, which are the the three- 
vector components of 



h kl L' kl Ss kl s' k $ Jkl j' kl (kki, hi, Ski\v Jkl \k' kl ,l' kl ,s' kl ). (31) 

and (• • • \v^ kl 1 • • • ) is the reduced relativistic two-body in- 
teraction. The two-body transition operators can be em- 
bedded in the three-particle Hilbert space: 



((kl)(m)\f k i(e k i)\(kl)(k)') := 



S(K, 



kl 



ki 



M, 



(0) 



(25) 



(26) 



where L(-) denotes a rotationless (canonical) Lorentz 
boost. These choices are appropriate for an "instant- 
form" representation. We use this representation for the 
purpose of illustration. Other choices give similar formu- 
las. 

The non-relativistic Jacobi moment are obtained if 
the Lorentz boosts L are replaced by the corresponding 
Galilean boosts. Both the relativistic and corresponding 
Jacobi momentum operators have identical spectra. 

For the three-body problem it is useful to successively 
pairwise couple the one-body Poincare irreducible repre- 
sentations to obtain three-body Poincare irreducible ba- 
sis vectors. This is done using pairs of Poincare Clebsch- 
Gordan coefficients. There are three bases that differ in 
the order of the pairwise coupling. We use the following 
shorthand notation for the basis that couples the irre- 
ducible representation associated with the pair (kl) to 
the representation associated with the "spectator" parti- 
cle, m: 

\P,j,H, (kl)(m)) := \P,j,n;L k i,S k i,K k i,j k i,l k i,s k i,k k i). 

(27) 

Here K k i and k k i are magnitudes of the relativistic Jacobi 
momenta. They are related to the invariant masses of the 
two- and three-body irreducible representation by 



M 



(o) 

kl 



m, 



(28) 



h kl L' kl $s kl s> kl 8j M f M (kki , hi , Ski \t 33k (e k i)\k' k i , l' M , s' kl ). 

It is important to note that the interaction V k i is con- 
structed from vu rather that v k i ■ This is intentional; V k i 
commutes with Km rather than pk + pz- It cannot com- 
mute with both because the relation (|25|l between these 
variables involves kki , which does not commute with Vki ■ 
This requirement is used to establish three-body Poincare 
invariance. It violates cluster properties for the genera- 
tors, but cluster properties are still valid for the scatter- 
ing matrix |15| . 

The two-body interactions that appear in the Faddeev- 
Lovelace equations have a complex relation to the inter- 
action Vki introduced above. 

The first step in the derivation of the interaction terms 
in the Faddeev-Lovelace equations is to add the interac- 
tion Vjk to the (jk) subsystem in the presence of a third 



particle spectator, /. The mass operator M { 



(jk)(l) 



is de- 



fined by its matrix elements in the basis 127|) as 
<P',/, A1 ',(jA : )(0 , |M o - fc)(0 |P,j, Ai ,(jfc)(0) = 



S(P' - P)6 fj S, l ,((jk)(iy\M Ukm \(jk)(l)) (33) 



where 



((jk)(iy\M um) \(jk)(i)) 



K' jk K jk 



S L ' l 3 Js> s f Jj' 3]k (l' jk , s' jk ,k' jk \ 



and 



< )2 



+ Hi + 



m 



Kll- 



(29) 



\ljk, Sjk, kjk) 



(34) 
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where V^ k is the two-body interaction introduced in 
equation iTIUl) . 

It follows from the definition (|33[l that 



(0)1 



[M - fc)(0 ,P]_ =0, 



(35) 



and the matrix element of Mq-^)^) has no explicit P de- 
pendence in the basis (|27H . 

Two-body interactions Vu k )U) that appear in the 
Faddeev-Lovelace equations are defined by 



(36) 



These can be used to define a three-body mass operator 

M := M(°)+V- (12)(3) +V' (2 3 )(1) +V(3 1)(2) = (37) 

It follows from the definition (|33|l that 

[M,j(°>]_ = [M,P]_=0, (38) 

and that simultaneous eigenstates of M, P, y°> and j' ' -z 
transform as irreducible mass M spin eigenstates of 
the Poincare group. Since these eigenstates are complete 
on the three-body Hilbert space, a Poincare invariant 
dynamics is defined by the requirement that these eigen- 
states transform irreducibly with respect to the Poincare 
group. 

To construct dynamical equations for the three-body 
transition operators, first define residual interactions V a 
for each two-cluster partition a — (ij)(k) by 

V a = V-V a . (39) 
The three-body channel transition operators 

1 



T (W) := V + V a 



W — M + i0+ 



V" 



(40) 



where W is a constant, satisfy the Faddeev-Lovelace 
equations 



T ab (W) = V b + J2 Vc 



W-M c + i0+ 



—T C \W). (41) 



where the sum runs over two-cluster partitions, c = 
(jk)(l),(kl)(j) and (lj){k). The identity 



where 



1 



W - M c + iCH 



= f c {W) 



1 



W - M + i0+ 



f c (W) := V c + V c 



1 



W — M c + i0+ 



Vc 



(42) 



(43) 



leads to the equivalent form of the Faddeev-Lovelace 
equations: 



T ab {W) = V b + J2Tc(W) 



l 



W - Mo + i0+ 



T ch (W). (44) 



The difficulty with the relativistic Faddeev-Lovelace 
equations is that the transition operator T C (W) differs 
from the off-shell two-body transition operators. This 
can be remedied by the observing that the cigenfunctions 
of M12 and M( 12 )(3) differ by delta functions. Our goal is 
to replace these by expressions involving non-relativistic 
transition operators. 

It follows that on the right half shell, 



((jk)(iy\(M m{l) -M )\(jk)(iy- 



(45) 



where the ) ± indicates a scattering eigenstate with in- 
variant mass eigenvalue 



w = \/ K h + Ak % + 47712 + \l K h + m2 ( 46 ) 



and 



((my\T% k (e)\(m)) 



((jk)(iy\(M Jk -M Jk o)\(jk)(iy 



(47) 



for e = 2^Jk 2 k + m 2 . Since the initial states are eigen- 
states of M( jfc )(;) and M^ k respectively and the final 

states are eigenstates of M^ and M^) respectively, it 
follows that 



(ijmy\Ti jkm (w)\(jk)(i)) 



2jk% 



2-i /k' 2 k 



Am 2 



Am 2 



((my\TT(e)m)(i)) 



(48) 



for the right half shell transition matrix elements. There 
is a similar expression for the left-half-shell transition 
amplitudes. 

The right hand side of this equation can be expressed 
in term of the reduced half-shell two-body transition 
operators which are related to the corresponding non- 
relativistic transition operators by l|15fl i|17fl an d 122fl . 
Combining these equations with l|32(l gives 

((my\% km (w)\m(i)) = 



^L kl L' kl Ss kl S' kl Sj kl j> k[ 



KKu - K' u ) 



-kl^ k i 
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2jk 2 k +m 2 + 2Jk' 2 k +m 2 



K )k + 4k% + Am 2 + JKf k + Ak' 2 k + Am 2 



Am 



e(fc w ) +e(k' kl ) 
for the CPS method 



T^(KiJki^'ki\tiA k ki/m)\k kl ,l kh s kl ) (49) 



While the purpose of this paper is to discuss the conse- 
quences of avoiding this additional step, we discuss this 
step briefly for completeness. 

The first method for constructing the off-shell transi- 
tion operators from the half-shell operators uses the spec- 
tral decomposition which can be expressed compactly by 

(f\ft(W")\i) = 



((jWY\Tl jkm (w)\uW)) 



S(K kl -K' kl 

u L kl L kl <JS k iS kl Vj k ij kl K K i 



2jk 2 k +m 2 + 2Jk' 2 u + m 2 



"jfc 



K ]k + 4k % + 4 ™ 2 + v K % + Ak fk + 4 ™ 2 



— ( k 'ki > I'm > s 'ki \ tJ (e(kki))\ k k l,hl,Ski)\ —■ 
m V 171 

for the GLC method, and 

((my\% m (w)\um)) = 



5(K kl -K' kl ) 

VL hl L hl "S M S kl Uj u j k i K,,K' 
t\ k il\ kl 



(50) 



a[a 



W" -W b + i0+ 



(r i \Vt\a" l+ )(a"' + \Vi^ 
W" - W'" + i0+ 



(52) 



where \b) denotes the two body bound states (we as- 
sume here that there is only one corresponding to the 
deuteron). Subtracting (|52[1 from the same equation with 
W" -► Wi gives 

(f\ft(W")\i) = {f\ft(W i )\i)+ 



(f'\b){b\i)[W b -W f ,}[W b -Wi\x 



1 



W" -W b + i0+ Wi-W b + i0+' 



d[a"']</'|f^(W"')|a"')(a"'|^(W"')|i)x 



k)k 



kfk 



±k% 



Am 2 + 



Kfk 



±kf k 



Am 2 



[W" - W f .}[W"' - Wi]x 



1 



1 



h (k'ki)(<l(k' k i), l'kn s 'ki\ t3 nr{ e rA<l{k k i))\q(k k i), l ku s M )h{k k i) 

(51) 

for the GK method. Unfortunately these relations only 
hold for the half-shell transition matrix elements. The 
relations H49fl - I|50|) are exact. 

The fully off-shell matrix elements of T C (W) that are 
needed as input in the Faddeev Lovelace equation can be 
obtained from the two-body bound state wave functions 
and the half-shell transition matrix elements by quadra- 
ture, or directly from the half-shell transition matrix el- 
ements by solving the first resolvent equation. Either of 
these methods provides a means for constructing the ker- 
nel of the relativistic Faddeev-Lovelace equations from 
solution of the relativistic Lippmann-Schwinger equa- 
tions. This step, of constructing the fully off-shell in- 
put to the Faddeev Lovelace equation, is the most non- 
trivial difference between the non-relativistic and rela- 
tivistic three-body problems. It involves either an addi- 
tional spectral expansion or solving an additional integral 
equation. 



W" 



W" + i0+ Wi - W" + i0+ ' 



(53) 



which expresses the off-shell transition matrix elements in 
terms of the half shell matrix elements using the spectral 
decomposition. These methods were used in the original 
Glockle, Lee and Coester paper and by Glockle and Ka- 
mada [|J. We also used this method in our calculations 
of in the next section. 

An alternative is to use the first resolvent equation 
with the definition of T^ a (W) to shift the value of W 

(f>\fi°(W)\i) = 



Wf 



^w 2 +K 2 + ^w 2 +K 2 



Wf + w 



K 2 +y/w" 2 +K 2 



{f\ft(w f )\a")d[a"]x 
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Wf — W 



Uwj + Kl + ^Jw" 2 + iff) (W - W") 



(a"\f^(W)\i) 



(54) 



This uses left half shell two-body transition operators 
to calculate the off-shell input to the Faddeev-Lovelace 
equations. The left half-shell two-body transition matrix 
elements have similar relations to the non-relativistic left 
half shell two-body transition transition matrix elements 
as the right half off shell transition matrix elements. This 
method can be used to systematically study the size of 
the corrections due to the off shell effects. 

The kernel of the Faddeev Lovelace equations has the 
form 



((i / i , )(fe , )|f (i , i , )(fc/) (W)|(i'/)(fe')) 



W 



M> 



Equations (|49l) . (1501) or l|51|) give expressions for 
((i'j'W)\T{i>?)(k>)(W)\(iy')(k')), when W = W, in 
terms of the non-relativistic two-body transition oper- 
ator. 

The non-relativistic kernel has a similar structure. The 
recoupling coefficients ((ij)(k)\(i' j')(k')) , which relate 
different orders of coupling of Poincare irreducible repre- 
sentations are replaced by coefficients that relate differ- 
ent orders of coupling SU (2) irreducible representations. 
The matrix elements {{i' j')(k')\f {t , r){y) {W)\{i' j')(k')) 
are replaced by off-shell two body transition matrix el- 
ements, and the energy denominator — Mg) is re- 
placed by the corresponding non-relativistic energy de- 
nominator. 

The energy denominator and transition matrix ele- 
ments in the relativistic and non-relativistic case are re- 
lated by kinematic factors when T is on the right half 
shell; this is precisely the point where the energy denom- 
inator has an integrable singularity. This suggests that 
the off shell contributions for T might be suppressed in 
the dynamical equations. 

Relativistic effects also appear in the permutation op- 
erators ((ij)(k)\(i'j')(k')}. The difference between the 
relativistic and non-relativistic coefficients can be easily 
seen when the recoupling coefficients are computed using 
the Balian-Brezin method j2f| [13 • 

The non-relativistic overlap matrix that relates the 
(12)(3) coupling scheme to the (23)(1) coupling scheme 
is: 

(P,j,/x ! (23)(l)|P'iV(12)(3)) = 



S(P-P')6(E-E 



2j + 1 



W=-i 



v ^ 87T 2 (mi + to 2 )(toi + m 3 ) 



Y£ 2amL (K 23 )y z ; 3mi (k a3 )y iiami (K 12 )r, i2 „ i; (k 12 ) x 



(j, /i"|L 2 3, m Ll S23, ms) (S23, m s \si, mi, j 23 , ^23} X 



O23, ^23^23, mi,S23, m 23 |s 2 , m 2 , s 3 , m 3 ) x 



(Sumu\si,mi,82, m 2 )(ji 2 ni 2 \l' 12 , m 'u si2,m 12 ) x 



(S 12 m s \s3,m 3 , ji2, Vi2){J, (J- \L 12 ,m L ,S 12 ,m s ) (55) 

where the angular variables can be computed for any 
convenient set of vectors of ki 2 , K i2 , k 23 , K 23 that have 
the correct values of k 2 2 , k 23 , K\ 2 and K\ 2 . 

An identical computation can be done in the relativis- 
tic case with the following result: 

(P,i,/i ) (23)(l)|P'iV(12)(3)) = 



S(P -P')S(M^ - M< »')|^ J2 £ 



M "=—3 



8tt 2 



^12^23^12-^23 



(Wl(/Cl 2 ) + U} 2 (k 12 )Y 



Wl (fel2 ) W2 (fcl2 )u>3 (if 12 ) Wl2 (if 12 ) 



1/2 



(uj 2 (k 23 ) + W 3 (fc 23 )) 3 



-,1/2 



U) 2 (k 23 )u) 3 {k 23 )uil (if 23)w23(if23) 



Y L 23 m L (K 23 )y i * 3 „ li (k 23 )Y L;2m ^(Ki 2 )y' J / i2m /(ki 2 )x 



(J,(i"\L,m L , S, m s )(S, m s |si, mi, j 23 , /x 23 ) x 



O23 1 M23I h mi, s 23 , to 23 ) (523, rn 2 3 \ s 2 , m 2 , s 3 , m 3 ) x 



^xrftx [R»(L(Jt u ), k' 12 )]^ 3 ^ 3 [i?- 1 (L(K 23 ), k 32 )] x 



(5i2mi2|si, mi, s 2 , m 2 )(jx 2 ni 2 \l', m[, si 2) m 12 ) x 



(S'm' s \s3, m 3 , j'12, H12) (J, fJ."\L', m' L , S' , m' s ) (56) 

Comparison of the two coefficients shows three differ- 
ences; 

1. The energy conserving delta function 6(E — E') in 
(|55H is replaced by a delta function in the invariant 
mass a(M<°> - Ml )'). 
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2. The second modification is that in the relativistic 



expression the coefficient 

87r 2 (mi + m 2 )(m 2 + m 3 ) 



TO 1 m 2 m 3 



(57) 



is replaced by 



8vr 2 



k 12 k23K 12 K: 



2.; 



(wi(fci 2 ) + w 2 (fci 2 )) 3 



Wi(fci 2 )a; 2 (A;i 2 )w3(i ; sri 2 )wi 2 (i ; iri 2 ) 



1/2 



(w 2 (fc 23 ) + w 3 (fc 23 )) 3 



-,1/2 



^2(^23)^3(^3)^1(^23)^23(^23) 



(58) 



which agrees with the non-relativistic expression in 
the limit that all of the relative momenta vanish. 

3. The third modification to the recoupling co- 
efficients involves the addition of momentum- 
dependent rotations that appear between the 
SU(2) Clebsch-Gordan coefficients. 

This choice of recoupling coefficient assumes that the spin 
in the relativistic problem is the canonical spin. In the 
relativistic expression R w (A, p) is the Wigner rotation 



R w (A,p) :=L- 1 (A P )AL(p) 



(59) 




where L(p) is the rotationless Lorentz transformation 
boosts a particle at rest to momentum p. The subscript 
K 23 indicates that the starting four vector is 

(w 23 ,K 23 )) = 



(60) 

As with the multiplicative factors, when all of the mo- 
menta vanish these rotations become the identity, and 
the expression reduces to the non-relativistic expression. 

In the Balian-Brezin form the transition from the non- 
relativistic to relativistic recoupling coefficients involves 
including some simple momentum dependent factors. 

The form of the relativistic recoupling coefficients will 
depend on the specific representation of the relativis- 
tic dynamics. The coefficients given above are in an 
"instant-form" representation. 



V. RESULTS AND DISCUSSION 

In this section consider the contributions of various 
terms to the kernel of the relativistic Faddeev Lovelace 
equations. 



The most complex modification of the non-relativistic 
Faddeev Lovelace equations involves the treatment of the 
off-energy shell kernel. As was mentioned in the previous 
section, in order to compute this from a non-relativistic 
two body transition operator fit to data, it was neces- 
sary either to utilize a full spectral expansion or to solve 
an integral equation. Both of these procedures involve a 
significant effort beyond that required to solve the non- 
relativistic equation. On the other hand we showed that 
the half-on-shell kernel could be calculated exactly in 
terms of the non-relativistic transition matrix elements 
by introducing only multiplicative kinematic factors. In 
addition, the on shell value was closest to scattering sin- 
gularity in the matrix elements of the three body free 
resolvent. 

The first test that we perform is to investigate the im- 
pact of simply including the kinematic factors that are 
needed to construct the half-on-shell kernel. These ap- 
proximations are given by equations I|50|). and IpTTJl. 

To set these approximations we consider a simple 
model without spin. The non-relativistic dynamics is 
given by the Malfliet-Tjon V interaction which has 
a long range attractive part and a short range repulsive 
core. Using this model we compare 

1. The exact calculation of off shell matrix elements of 
T(W) computed using the spectral expansion lf53|l . 

2. These approximations of T(W) given by equations 

(|50"jl . and lj5l]) . and by keeping only the first 
term in Eq. Ij53|l . ignoring the dispersive contribu- 
tions. 

3. The non-relativistic off shell transition operator. 

The following kinematic variables were varied in these 
comparisons: 

1. W = 3m + Enn (three-body energy); 

2. Kjj (spectator momentum); 

3. k (initial two-body momentum); 

4. k' = k (initial=final two-body momentum). 

We also considered all three methods of using the non- 
relativistic two-body dynamics as input. 

The results of these calculations are shown in figures 
3 - 14. The plots compare the difference between the 
non-relativistic and exact results divided by the exact 
result. This is done for each method. Figures 3-6 use 
the CPS method, figures 7-10 plot the same quantities 
for the GLC method, and figures 11-14 plot the same 
quantities for the GK method. The first plot in each set 
varies the off shell energy for initial and final values of kij 
equal and for a fixed value of Kij . The next plot fixes the 
energy Kij and varies the initial and final values of fey . 
The third plot varies the final kij keeping everything else 
fixed and the fourth plot varies Kij holding everything 
else constant. 
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FIG. 3: Fractional deviation of approximate T matrix com- 
pared to exact evaluation of Eq. 1531 using the CPS method, 
as a function of the 2-body energy (MeV) in the 3-body CM 
frame. Dotted line uses approximate value described in the 
text. Dashed line uses the nonrelativistic value. 

In the CPS and GLC cases the error in the approxima- 
tions l|49|) or H5U|) to the exact off shell matrix elements is 
no more that a few tenths of a percent in all four cases. 
They are significant improvements over the straight non- 
relativistic result. In the GK case the errors can be as 
much as 5% percent in some kinematic regions, although 
the errors are typically about twice the size of the errors 
using the other two methods. 

These calculation suggest that the bulk of the rela- 
tivistic effects in the transition operators in the Faddeev 
Lovelace kernel comes from the kinematic factors that 
appear in (J201 or (|5l|) 

The presence of the free resolvent in the Faddeev 
Lovelace kernel only serves to improve the approxima- 
tion to the full kernel, since it is largest on the half shell. 

The other factor that appears in a full relativistic three 
body calculation are the three body recoupling coeffi- 
cients. These differ from the non-relativistic coefficients 
by kinematic factors and the presence of momentum de- 
pendent spin rotation functions. In many relativistic cal- 
culations the overall contribution of the spin effects have 
proved to be unimportant, however most of these appli- 
cations have involved bound states. 

It is a simple matter to compute relativistic and non- 
relativistic recoupling coefficients. To compare them the, 
delta function in either invariant mass or non-relativistic 
internal energy is replaced by a delta function in the mo- 
mentum followed by the appropriate Jacobian factors. 
The comparison can then be made between the coeffi- 
cients of the delta functions in the two cases. Figure IVTll 
is typical of the behavior of the recoupling coefficient. 
The curves are computed for k\2 = K23 = K\i for a 



FIG. 4: Fractional deviation of approximate T matrix com- 
pared to exact evaluation of Eq. 153H using the CPS method, 
as a function of the 2-body relative momentum (fm -1 ), initial 
and final momenta kept equal. Dotted line uses approximate 
value described in the text. Dashed line uses the nonrelativis- 
tic value. 



representative set of spin quantum numbers. The solid 
curve is the ratio of the non-relativistic to the relativistic 
recoupling coefficient when they are made to have the 
same delta function. The dotted curve shows the ratio of 
relativistic coefficient with the rotation functions turned 
off to the full relativistic recoupling coefficient. 

The curves show that the non-relativistic curve is 
about 2% higher than the relativistic curve at 0.5 GeV. 
It grows to just over 5% at 1 GeV. Including the correct 
multiplicative kinematic factors without including spin 
rotations leads to considerable improvement at 0.5 GeV; 
and it grows to just over 2% at 1 GeV. This is consistent 
with previous calculations that suggest the the recoupling 
coefficients have small effects. 

With the Balian-Brezin method the inclusion of both 
the spin factors and kinematic factors involve minor mod- 
ifications to an existing non-relativistic program. 

Our conclusion is that negligible errors will be made if 
the relativistic Faddeev-Lovelace kernel is approximated 
by replacing the transition matrix elements of T by the 
expressions in (|49|l . (|50|) . and (J5TJ- Care must be used 
to match the method used to the fitting procedure used 
to determine the non-relativistic two body interaction. 
The recoupling coefficients arc not needed for calcula- 
tions dominated by lower momenta; and they lead 5% 
effects at 1 GeV. The inclusion of the correct kinematic 
factors reduces errors considerably. 
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FIG. 5: Fractional deviation of approximate T matrix com- 
pared to exact evaluation of Eq. 1531 using the CPS method, 
as a function of the 2-body initial relative momentum (fm~ ). 
Dotted line uses approximate value described in the text. 
Dashed line uses the nonrelativistic value. 



FIG. 7: Same as Fig. 3 except using the GLC method. 
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FIG. 8: Same as Fig. 4 except using the GLC method. 



FIG. 6: Fractional deviation of approximate T matrix com- 
pared to exact evaluation of Eq. 1531 using the CPS method, 
as a function of the spectator momentum (fm _1 ). Dotted 
line uses approximate value described in the text. Dashed 
line uses the nonrelativistic value. 



VI. COMPARISON TO APPROACHES BASED 
ON QUASIPOTENTIAL THEORY 

The development and the results discussed above uti- 
lize Wigner's formulation of relativistic quantum me- 
chanics, using a direct construction of a unitary rep- 



resentation of the Poincare group on the two-nucleon 
Hilbert space . The two-nucleon interactions have spe- 
cific connections both to the Schrodinger equation and 
to the three-nucleon problem. One can also approach 
these problems within the framework of a quasipoten- 
tial theory. The first results for three-nucleon binding 
were obtained by Stadler and Gross Il2|. using the so- 
called spectator approximation (23, |2J, |2|| . Within this 
framework they can obtain the observed triton binding 
while maintaining a fit to the observed nucleon-nucleon 
phase shifts and the deuteron binding. When they com- 
pare the results of their "full" calculation to a nonrela- 
tivistic approximation, they find that most of the new 



13 





GLC Kinematics; k f = 1.0; E = - 


10.0; P = 1.0 


0.005 










A(NR"Efac- full l/full 

A(NR-full)/full 




0.0 












H 






■0.01 






-0.015 






-0.02 








0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 

k, 



GK Kinematics; k r = k, = 1.0; P = 1.0 



-0.U02 
■0.U04 



■0.008 

-0.01 



FIG. 9: Same as Fig. 5 except using the GLC method. 



FIG. 11: Same as Fig. 3 except using the GK method. 
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FIG. 10: Same as Fig. 6 except using the GLC method. 



FIG. 12: Same as Fig. 4 except using the GK method. 



effects come from the inclusion of negative-energy states. 
These would correspond approximately to three-nucleon 
interactions in Poincare invariant quantum theory. They 
also find that so-called "boost effects" in their framework 
are small and give a repulsive contribution to the triton 
binding. Those results are consistent with the results 
obtained in Poincare invariant quantum theory. 

An earlier work by Sammarruca and Machlcidt 13] 
examined the effect of kinematic factors within frame- 
works based on field theories. They made use of an ear- 
lier relation obtained by Brown, Jackson and Kuo |2t| . 
who used quasipotential methods to argue that the non- 
relativistic interaction should include non- localities. The 



Brown, Jackson, Kuo method leads to a modified inter- 
action, however unlike the methods discussed above the 
new interaction must be refit to the phase shift data. 
This method leads to the following relation between their 
minimal relativistic and nonrelativistic t matrices: 



m V ni 



(61) 



Sammarruca and Machleidt then use this modified two- 
nucleon interaction within their relativistic framework 
and get a net attractive correction to 3 H binding when 
using this factor. Again, this is understandable at least 
naively when considering that the t matrices receive es- 
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FIG. 13: Same as Fig. 5 except using the GK method. 
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FIG. 14: Same as Fig. 6 except using the GK method. 



sentially the opposite factor from those discussed above 
[Eqs.E 



3], However, both the Brown, Jackson and Kuo 
method and the CGL method in principle require a refit- 
ting of the phase shifts to get an exact phase equivalence. 
Thus we cannot make a direct comparison to these meth- 
ods. 



VII. CONCLUSION 

We have explored here the quantitative impact of the 
requirement of relativistic invariance with the framework 



FIG. 15: Ratio of the recoupling coefficients. 



of models discussed here, the dominant effects are mani- 
fested in kinematic terms that can easily be incorporated 
into a Schrodinger-based momentum-space calculation, 
namely, multiplicative factors of the E/m form, and en- 
ergy denominators that employ the relativistic energy- 
momentum relation. There are additional small effects 
that come from relativistic corrections to the recoupling 
coefficients. These can be minimized by including the 
correct relativistic kinematic factors. A correct treat- 
ment of these terms involves only minimal modifications 
of the non-relativistic recoupling coefficients. Quantita- 
tively the effects of all of these corrections are small, and 
contribute repulsive corrections to three-nucleon binding. 
This means that the resolution of the discrepancy be- 
tween Poincare invariant three-nucleon calculations and 
the experimentally observed binding energies must come 
from three-nucleon forces. The quantitative relativis- 
tic effects grow with internal momenta; non-relativistic 
three-nucleon scattering calculations at energies of sev- 
eral hundred MeV may well require corrections beyond 
the simple square-root factors. 
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